Screening of the Raman response in multiband superconductors - application to iron-pnictides 
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We performed model calculations of Raman responses for multiband 2D superconductors. The multiband 
effects of screening in the A\ g symmetry channel were investigated analytically and numerically for a band 
structure model mimicing ARPES data on iron-pnictide materials. An acceptable agreement between our model 
calculations and recent experimental data is demonstrated by modification of the band structure parameters. 
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I. INTRODUCTION 

A large variety of compounds belonging to the iron-based 
superconductors have been found since the discovery of super- 
conductivity in LaFeAsO\- x F x in the beginning of 2008 1 . 
They can be classified into four different families, divided by 
stoichiometric differences, exhibiting superconducting (SC) 
behavior induced by hole-doping, electron-doping as well as 
solely the application of high pressure (for a review see^. All 
compounds are layered materials with a common layer of an 
iron-pnictide in between different spacing layers or simply 
layers of FeSe. So the location of the SC inducing physics 
is the tetragonal Fe-containing layer leading to highly two di- 
mensional (2D) electronic structure 3 . Angle resolved photo- 
emission spectroscopy (ARPES) finds several bands crossing 
the Fermi surface 4 6 (FS) which is a novelty since, with the 
exception of MgB2, superconductors have been most often 
treated as single band materials. The mechanism that leads 
to superconductivity including the pairing symmetry of the 
SC order parameter is still under debate (a theoretical review 
is provided by RefP). Gap symmetries both with nodes (p- 
wave or d-wave) and without nodes (s-wave) are still under 
consideration^. 

Electronic Raman scattering (for a review see^) could clar- 
ify some of those unresolved issues. In the light scattering 
experiment different parts of the Brillouin zone (BZ) can be 
emphasized with different weights through the selection of 
incident and outgoing light polarizations with respect of the 
crystallographic directions. Thus by manipulating the direc- 
tions of the polarization vectors of the incident and scattered 
photons in a Raman experiment, it is possible to determine 
the symmetry of the pair-breaking excitations related to the 
symmetry of the SC order parameter 8 . Furthermore, the low 
energy tail of the Raman responses can be used to distinguish 
between a fully gaped FS or a SC gap with nodes allowing the 
excitation of quasi particles at arbitrarily small energies. 

The theory of electronic Raman spectroscopy in 
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superconductors 9 has be en v ery successful in applica- 
tions to high-T c cuprates^ED and now allows to make 
predictions for the Raman responses of the iron-pnictides. 
Recently model calculations of expected Raman responses 
for a superconductor with multi-band order parameter and 
including the long-range Coulomb interactions responsible 
for the multi-band screening effects have been performed 
by Boyd et. alP^. These calculations used multiple circular 
FS sheets as a band structure and for the majority of the 
work the Raman vertices are obtained by an expansion into 
FS harmonics for cylindrical FS. Taking the screening of 
the charged background into account but leaving out vertex 
corrections capturing effects of final state interactions a 
variety of different gap distributions on the FS sheets with 
and without nodes were used. Thus many different forms 
were obtained for the Raman responses which show the 
possible scenarios expected from the data. By comparing 
the singularities in different symmetry channels of such 
calculated Raman responses with the peak structure of actual 
Raman data the reconstruction of the momentum and fre- 
quency dependence of the superconducting order parameter, 
critical to understanding of the superconducting mechanism 
in these materials, is possible. In another stud y 12 ! 13 ! it has 
been established by analyses of the vertex corrections that 
for an A\ g extended s-wave gap symmetry, the A\ g Raman 
response has a true collective mode resonance peak below the 
fundamental gap 2Ao- Since the collective mode is predicted 
only for this particular pairing symmetry, conscientious polar- 
ization resolved low-frequency Raman scattering study of the 
iron-pnictide materials may provide a way to unambiguously 
distinguish between various suggested gap symmetries. 

Very recently the first electronic Raman scattering experi- 
ments on Ba(Fei- x Co x )2As2, a member of the 122-family, 
have been performed 14 showing a change of the Raman re- 
sponse in the SC state only in the A lg and B 2g channel. A 
clear peak of fog-singularity shape in the latter one and a 
broader peak in the former one at a slightly higher energy are 
observed and the low energy shape is interpreted as evidence 
for accidental gap nodes. 

In this paper model calculations of Raman responses for a 
multiband free electron model are performed. In chapter II 
the general theory of Raman scattering in the SC state is pre- 
sented and explicit forms of the Raman vertices and an ana- 
lytical expression for the Raman susceptibility are calculated 
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for a 2D tetragonal symmetry and a constant SC gap. Then 
the screening correction is taken into account adding an extra 
term in the A\„ channel. Multiband effects in this symmetry 
channel for two and more than two bands with an example of 
a four band response are presented in chapter III. The general 
effects found in chapter III are used in chapter IV to explain 
the Raman responses for a band structure model inspired by 
ARPES data. First a constant gap and then a /c-dependent gap 
of extended s-wave symmetry is used for which the single 
band structure is investigated with an only angular dependent 
model gap. Finally the Raman response with the extended s- 
wave gap is compared to experimental data and it is shown 
that the band structure model can be modified to come to ac- 
ceptable agreement. It is important to note that final state in- 
teractions leading to vertex corrections are neglected in this 
work. The consequences of including vertex corrections that 
could give rise to resonance modes or notably modify the Ra- 
man response are discussed in Refsj 12 l 15 l However, this work 
focuses on multiband screening effects and features of vertex 
corrections are left for future publications. 

The fact that in the applied assumptions analytical calcula- 
tions are possible allows us to find the origin of the obtained 
general mutliband effects in the screened A% g channel. With 
that the shape of the numerically calculated Raman responses 
for the particular band structure with extended s-wave symme- 
try can be explained and further applications to other multi- 
band SC might follow. Furthermore the nature of the main 
features of the single band Raman responses with the extended 
s-wave gap are identified analytically with a simplified model 
gap. Thus our work adds further details to the previous study 
by Boyd et. alP. 



II. THEORY OF ELECTRONIC RAMAN SCATTERING IN 
THE SUPERCONDUCTING STATE 

The Raman process is inelastic scattering of light in which 
the energy shift u> = ujj—ujg between the incident photon (luj, 
kf) and the scattered photon (ug, kg) is measured. In the case 
of electronic Raman scattering this two photon process creates 
charge density fluctuations in the sample and trough a virtual 
interband transition (in the non resonant case) an electron is 
excited from a state below the Fermi energy (Ep) to a state 
above. In a clean system the energy difference between those 
two states corresponds to the energy loss of the photon. It can 
be shown that the cross section of the scattered light is pro- 
portional to the imaginary part of the Raman susceptibility 10 
which is defined in linear response theory in the following way 



XRaman{q,t) = T Tr 
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[PqW,P-q(0)] (1) 



Here q = kj — kg, (3 = kgT (kg is Boltzmann's constant), 
Z is the partition function and Hq is the Hamiltonian of the 
unperturbated system. The perturbation Hamiltonian H' is 
proportional to the effective density operato^ 
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which contains fermionic annihilation (c n k ) and creation 
( c « k) operators in the n th band. Furthermore p q is propor- 
tional to the so called Raman vertex j n (k) which depends on 
the polarization direction of the incident and scattered pho- 
ton and thus on the experimental setup (light polarizers). This 
quantity will be specified and discussed in detail in the next 
paragraph. With the wavelength of the incident and scattered 
photon being in or near the visible spectrum the correspond- 
ing wave vectors and therefore also their difference q is much 
smaller than ^ A (2 A is the SC gap and vp q is the compo- 
nent of the Fermi velocity along q). This makes the q — » 
limit a reasonable simplification that leads to the explicit form 
of the Raman susceptibility in the SC state 

^nCl^ '") = JfzZ^n( k )^n{k,^)- (3) 
k 

The quantity A„ (k, uS) is the Tsuneto function^ 
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with the SC gap A„ (k), the SC quasi particle energy E n (k) = 
\J e^(k) + A^(k) and a small broadening factor a. This func- 
tion represents the creation of quasi particles out of the SC 
condensate. The first term is proportional to a SC coherence 
factor and is unity at the Fermi momentum. The further away 
from Ep the quasi particles are created the smaller this factor 
gets so the major contribution to the Raman response origi- 
nates from the vicinity of the Ep. The second term is a sta- 
tistical factor that takes into account already thermally excited 
quasi particles. Since the following calculations are made un- 
der the assumption that T = this factor is simply unity. The 
third term takes care of the necessary energy conservation and 
excludes a creation of quasi particles inside the SC gap. 

Now the Raman vertex will be specified by an analytical ex- 
pression. In the nonresonant case, where the electronic excita- 
tion does not include a real intermediate state, the Raman ver- 
tex can be calculated with an effective mass approximation 10 
in the following form 
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with to denoting the electrons mass and i and e 7 a being 
the light polarization vectors of the scattered and incident pho- 
ton, respectively. In this approximation the Raman vertex is 
proportional to the second derivative (curvature) of the band 
energy dispersion in the particular fc-direction selected by the 
two light polarization vectors which is proportional to the in- 
verse effective mass of the band in k ■ p - theory. Taking into 
account that a periodic crystal transforming according to a 
special point symmetry group is investigated group theoret- 
ical arguments can be used to further simplify the expression 
in Eq. j^. Any operator of the form Mj A — M^efe s 
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(i=initial, /=final) can be decomposed according to the point 
symmetry group of the crystal which will be taken to be D^h 
in this case. This leads to the decomposition of Mfi into 
projected operators O v corresponding to the irreducible rep- 
resentation v of the point groujP. 



M 



1 



1 



2 Alg (ef e% + 44) + -0 Blg (ef e% - 44) 



(6) 

which results, together with Eq. in the Raman vertices for 
the four symmetry channels 
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The different symmetry channels can be investigated by 
choosing certain combinations of the light polarization vec- 
tors through the setup of light polarizers in the experiment. 
Using an effective mass approximation (Eq. |5]l) to obtain the 
Raman vertices of Eq. (|7]) is a different approach than the 
one in Boyd et al. 11 where an expansion into FS harmonics is 
done leading to angular dependent vertices. In this work a 2D 
free electron model for the band structure will be used which 
results in a constant vertex for the Ai g channel since the sec- 
ond derivative of the energy band dispersion is proportional 
to the constant inverse effective mass -\ of the band. For 
a parabolic band all other symmetry channels are vanishing 
which is a consequence of the effective mass approximation 
in combination with a free electron model. This leads to the 
focus of the following investigations being on the fully sym- 
metric A lg channel. A tight binding model, which is the next 
possible step, will not lead to vanishing vertices for the non 
fully symmetric channels (except A^ g , the symmetry of the 
antisymmetric tensor^, which is usually negligible). 

Using a 2D free electron model for the energy dispersion of 
the n th band, which is determined by the effective mass m* 
and the chemical potential fj, n , the Raman susceptibility x-yy,n 
can be calculated analytically. First the Tsuneto function of 
Eq. Q is inserted into Eq. Q, the limit of a —> is taken and 
the sum is replaced by a 2D fc -integration. Evaluating only the 
imaginary part of the integral and inserting the result into the 
Kramers-Kronig relation leads to the analytical expression 

X77 ,„H = N F , n ( 7 „) 2 arCS ^li = N F , n { ln f F n (8) 



with the definition 
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and with x n = and Nf,ti being the 2D density of states 
(DOS) at Ep which is a constant proportional to |m* | (de- 
tailed calculation can be found irM^I). Since the cross section 
of the Raman process or Raman efficiency is proportional to 
7ro[x77,n] the one band response is obtained by taking the 
imaginary part of Eq. ([8j. For < x n < 1 (0 < cj < 2A n ) 
/m[arcsin(cc„)] is vanishing and the square root is entirely 
real leading to a vanishing Raman response. This is due to the 
fact that no quasi particles can be created without breaking a 
Cooper pair which has a binding energy of 2A n . At oj = 2A„ 
{x n = 1) the denominator hits zero and the function diverges 
as a square root singularity, a consequence of the diverging SC 
quasi particle DOS at this energy. Above u — 2A„ (x n > 1) 
the square root is entirely imaginary and in im[x 77i „] only 
i?e[arcsin(x n )] remains which is a constant of s leading to 
7m[x 77 .„] = f i?e[ — ^ 2 J for the single band Raman re- 
sponse. 

So far it was not taken into account that in the solid a 
charged background of electrons is present that reacts through 
long range Coulomb interaction on the charge density fluc- 
tuations imposed by the Raman process. The electron back- 
ground rearranges in a way to screen the perturbative elec- 
tric fields and thus works against the Raman process. This 
screening effect can be included in the theory within a ran- 
dom phase approximation like sum 10 that leads in the limit 
of small q (compared to the inverse coherence length £ and 
the Fermi wave vector kp) to an additional term in the Ra- 
man susceptibility. This term includes a k sum over the 1st 
BZ of 7(k) • A(k) that vanishes if 7(k) is not fully symmet- 
ric hence exhibiting sign changes in between different parts 
of the 1st BZ leading to those parts canceling out each other. 
Thus only the fully symmetric A\ g channel gets screened and 
the screened Raman susceptibility is 
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The second term in Eq. ( 10 1 is called screening term and con- 
sists of a normalization by a sum over terms shown in Eq. 
( p"2| ) and the square of a sum over terms shown in Eq. ( [TT) . 
Considering only one free electron band the Raman response 
vanishes. 
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This means that the screening of a free electron background 
is perfect and all charge density fluctuations are completely 
compensated by the rearrangement of the free electrons. 
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in. MULTIBAND EFFECTS IN THE SCREENED A lg 
CHANNEL WITH CONSTANT SC GAPS 



A. Two bands in A\ g symmetry 

In contrast to the rather trivial one band case interesting fea- 
tures appear as soon as two free electron bands are considered. 



Taking a closer look at the second term in Eq. ( 10 1 reveals 
that through the square in the numerator cross terms of the 
complex functions x-yi,n containing real and imaginary parts 
corresponding to distinct bands contribute to the two band Ra- 
man response. Having two bands furthermore adds complex- 
ity by the possibility of choosing different effective masses 
(magnitude and sign) and SC gaps. The general form of the 
screened A\ g Raman response shown in Eq. (p~3l) extended for 



two bands can be transformed into a simpler form 

Im [XRamani^)} = Im [iV>, 171-^1 + iYp^f-Fb] 
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For the sake of generality the effective masses (determining 
7„ and N Ft7l ) and the SC gap (present in F n ) are taken to 
be different for both bands performing the transformation in 
Eq. 



(14i 



We note here that for a two band model the vertex 
corrections due to final state interactions can notably modify 
the Raman response (compare Eq. ( 14 1 to Eq. (6) in ReffSJ). 
The first thing that becomes obvious in the form of the last 
line of Eq. ( 14 1 is that for equal effective masses (equal 7„'s) 
the two band Raman response is vanishing independent of the 
choice for the two SC gaps. So the screening is again perfect. 

Further investigations of the different two band features re- 
quire a case by case study of all possible combinations of the 
two tunable quantities, namely the effective mass and the SC 
gap. First the case of two equal SC gaps (Fi — F 2 — F) is 
considered. In order to see the effect of screening in this case 
one has to take a closer look at the screening term which is 
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with 
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M* 



(16) 



In Eq. (15 1 all common constants are absorbed in C. The 



absolute values of the effective masses in Eq. ( 16 1 have their 
origin in the fact that the DOS is always a positive quantity in 
contrast to the Raman vertices for which the sign of the effec- 
tive mass remains. In the numerator of the inverse effective 
screening mass the absolute value of the effective masses 
cancels out leaving only a sign dependence. This feature has 
its origin in the two dimensionality of the band structure for 
which the DOS is a constant proportional to |m* | and will not 
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FIG. 1: Comparison of unscreened (black) and screened (red) Raman 
responses in Ai g symmetry for two free electron bands with unequal 
Raman vertices. In the left column the case of two equal constant 
SC gaps A„ and opposite signs of the effective masses is shown for 
which the screening is vanishing. In the right column the case of 
two different SC gaps is displayed with different combinations of 
effective masses mj and m\ for the two bands. In all three cases the 
singularities at twice the SC gap values are removed. Color online. 



remain for ID or 3D systems where the DOS includes an en- 
ergy dependence. For two bands with equal signs this leads 
to a finite screening term and thus a partial screening (if the 
effective masses have different absolute values) while for un- 

in Eq. ( 16 1 is zero. Thus for 



1 

A/* 



equal signs the numerator of 
a combination of a hole-like and an electron-like band with 
equal SC gaps the screening term vanishes and the two band 
Raman response is entirely unscreened no matter how differ- 
ent the absolute values of the effective masses of the two bands 
are. This case is displayed in the left column of Fig. [T] where 
the unscreened Raman response (black) and the screened one 
(red) show no difference demonstrating the vanishing screen- 
ing term. Additionally it has to be noted that the square root 
singularity previously found in the unscreened one band Ra- 
man response remains in the screened two band Raman re- 
sponse for equal SC gaps. 

Now the case of unequal SC gaps is investigated. Reinsert- 
ing the explicit form of F n (see Eq. Q) into the simplified 
form of the screened two band Raman response (last line in 
Eq. 



14 1) results in 
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AS n — arcsin (x n ) 



(17) 



(18) 



The form in Eq. (17i shows that the singularities at to = 
2A12, present in the unscreened case, are both removed in 
the screened case since the square root terms SQi and SQ 2 
hit zero separately always leaving one of them finite. This 
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case is illustrated in the right column of Fig. [T]with differ- 
ent combinations of effective masses for the screened Raman 
response (red, orange, yellow). For equal effective masses (or- 
ange) a broad dome, peaked in the middle of twice the two SC 
gap values, appears. When both bands have different effective 
masses (red, yellow) this shape changes into a sharper peak 
located in the vicinity of twice the SC gap corresponding to 
the band with the smaller effective mass. 



B. More than two bands in Ai g symmetry 

Having seen the interesting features originating from the in- 
terplay of two bands through their cross terms in the screening 
term the next logical step is to investigate what happens by in- 
cluding more bands into the Raman response. Writing down a 
general form for more than two bands (analog to the upper two 
lines of Eq. ( fT4) i) it becomes clear that the screening vanishes 
only in very special cases. Only if all SC gaps are equal and 
an even number of bands consisting of the same amount of 
hole-like and electron-like bands is considered the Raman re- 
sponse will be unscreened independent of the absolute values 
of the effective masses. In all other cases the strength of the 
screening will depend on the combination of effective masses 
of the contributing bands. A transformation analog to the one 



done in Eq. ( 14 1 including n free electron bands leads to 



Im [XRaman(u)} = Im 
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Here i,j and / run from 1 to n. This simplified form shows 
that the screened n band Raman response consists of a sum 
over all possible combinations of two band terms (with i ^ j) 
each normalized by all n bands. Furthermore it is obvious that 
terms including two bands with equal Raman vertices vanish 
and that the dominating terms in the sum are the ones includ- 
ing two bands with opposite signs since for them the Raman 
vertices in the bracket add up. 

Now the issue of the removal of the singularities previously 
seen for the two band response with unequal SC gaps will be 
studied for more than two bands. In order to simplify this 
investigation the example of a four band Raman response, in- 



stead of the general n band form in Eq. ( 19 1, will be consid- 
ered first. For four bands there are six possible combinations 
of two band terms that add up in the sum of the total Raman 
response. Each of those six terms can be written as 



Im 



with 
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If the two bands i and j in Eq. ( 20 1 have different SC gaps 
SQi and SQj will hit zero at different energies. Thus the term 



Li = SQj will remain finite when Lj = SQi = and vice 
versa leading to a finite denominator of Eq. ( |20] > for all ener- 
gies. For equal SC gaps in Eq. (20i SQi and SQj are equal 



and will vanish at the same energy resulting in a vanishing de- 
nominator and hence a singularity at that energy. This means 
that as soon as there are two bands with equal SC gaps con- 
tributing to the four band Raman response there is one term in 
the sum that is singular at twice this SC gap. An illustration of 
this is seen in Fig. [2]where different cases for the SC gaps are 
considered. The unscreened Raman response (black) shows 
the location of the singularities and the screened one (red) dis- 
plays which ones get removed through screening. Whenever a 
gap is present more than once the screened four band Raman 
response contains a singularity because of the above given ex- 
planation. The generalization of the four to the n band case 
can easily be done by adding further terms to the denominator 



of Eq. ( 20 1. The argumentation will remain the same as above 
and the discovered result that only singularities corresponding 
to uniquely present gap values are removed holds true. 



IV. APPLICATION TO A MULTIBAND FREE ELECTRON 
BAND STRUCTURE BASED ON ARPES DATA FOR AN 
IRON-PNICTIDE SUPERCONDUCTOR 

A. Band structure model 

Here the general investigations of free electron Raman 
responses in the SC state will be applied to a particular 
band structure model inspired by ARPES measurements. In 
Bai- x K x FeAs the band structure was obtainecP^consisting 
of two hole like, almost parabolic bands around the T-point 
and a propeller-like structure around the M-point. This pro- 
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FIG. 2: Unscreened (black) and screened (red) Raman response of 
four free electron bands with constant SC gaps A n in Ai g symme- 
try. Singularities corresponding to uniquely present gap values are 
removed by screening and the dome-like shape previously seen in 
the two band case appears. Color online. 



6 



peller has an electron-like parabolic band in its center which 
is surrounded by four hole-like blades. In order to mimic this 
structure within a free electron model the two bands around 
the r -point and the one around the M -point are approximated 
by a 2D parabola which has the form 



2m) 



((k x - g) 2 + (k y - h) 2 J - \ii 



(22) 



with g,h E Z. The blades are approximated by 2D ellipsoids 
having two different effective masses in the k x 
fc y -directions, represented by 
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is a direct consequence of the form of the Raman vertex and 
the 2D DOS, both different for parabolic and elliptical bands, 
is the vanishing screening for two bands of opposite signs and 
equal constant SC gaps. Since the Raman vertex of an el- 
liptical band in A\ g symmetry is not just proportional to the 
inverse effective mass but to the sum of the two different ef- 
fective masses in the particular fc-direction and the 2D DOS is 
proportional to the square root of the product of the two effec- 
tive masses the inverse effective screening mass differs from 
Eq. ( [To} . For two elliptical bands it changes to 



(26) 



M* 



(23) 



i&a, with & 1 e Z and & o = -1, 0, 1, 



Here 8 t = £^1 + ^7=^,2, wmi c,i,l t m aim £ i)2 

which is extracted from 5 . The effective masses and chemical 
potentials /i„ used for the different bands are given in TABLE 
IP. Eq. ( 22 1 and Eq. ( 23 1 with the parameters from TABLE 
Hlead to the band structure shown in the left of Fig. [3]and its 
corresponding FS displayed in the top middle of Fig.[3] 



where the plus describes two bands of the same sign and the 
minus two bands of opposite signs and m* x+y and m* x _ y for 
band n are denoted by m* 1 and m* 2 for a shorter notation. 
Eq. ( |26] > shows that for two elliptical bands with equal signs 
the screening is not vanishing (just like it is the case for two 
parabolic bands) and that its strength depends on the absolute 
values of the effective masses. Furthermore it can be seen that 
for two bands of opposite signs and different ratios between 
m* 1 and m* 2 the screening term remains finite while for the 
same ratio of the effective masses both terms in the numerator 



TABLE I: Effective masses and chemical potentials used for Fig. [3] of Eq. (|26]l cancel each other and the screening is vanishing. 



band 


inner T-par. 


outer T-par. 
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m* [m e ] 


-2.1 


-6 


0.8 


-0.61, -2.46 


ft [meV] 


35 


45 


-17 
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Combining Eq. (|7]i with Eq. ( |22] > and Eq. ( 23 1 results in 
explicit Raman vertices for the used band structure model. For 
the parabolic bands we obtain 



Ta 13 

7B lfl 



m 
'ml 

=7k 



TA 2a 



(24) 



and the Raman vertices for the elliptical bands are 





m 1 




= ^ \ 




m I 




= ^ \ 






-ft. 





1 



1 



"x+y 
1 

l x+y 



x-y 
1 

l x-y 



(25) 



For the used band structure all bands contribute to the fully 
symmetric A\ g channel while the Bi g channel only has a 
contribution from the elliptical blades and the other symme- 
try channels completely vanish. This is a consequence of the 
free electron model in combination with the effective mass 
approximation used inhere. 

At this point it is interesting what the new aspects are that 
come into play by including elliptical bands. The effect that 



B. Screened A\ g Raman response with a constant gap 

In order to be able to calculate the Raman response of the 
introduced band structure shown in Fig. [3] the distribution 
of the SC gap size along the different bands has to be spec- 
ified. As a first step the SC gaps are assumed to be constant 
on each band which makes it possible to directly apply the 
general investigations of chapter III. ARPES measurements 
find two different SC gaps 6 of which the smaller one is lo- 
cated on the outer T-parabola and the bigger one is present on 
all the other bands. The magnitude found for the smaller gap 
is roughly half the size of the bigger one. Thus the distribu- 
tion of the SC gaps used here consists of a smaller gap on the 
outer T-parabola and a bigger one on the inner T-parabola, 
the M-pocket and the blades both related by 2Ai = A2. The 
screened (red) and unscreened (black) A\ g Raman response 
for these assumptions presented in Fig. fflis reminiscent of 
the one shown in the lower right of Fig. [2] In both plots the 
singularity at twice the smaller SC gap is removed since it is 
uniquely present and the one at twice the bigger SC gap re- 
mains because three bands exhibit this same gap value. As 
previously mentioned the B 2g Raman response just has a con- 
tribution by the blades leading to a simple square root singu- 
larity structure like in the left hand side of Fig. [T] with the 
singularity located at twice the bigger SC gap. The other two 
symmetry channels have vanishing Raman vertices and hence 
show no response at all in this model. 
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FIG. 3: Left: 2D free electron band structure model based on ARPES datPl. Four different bands are crossing the FS (red plane). The inner 
r parabola (blue) and the outer T parabola (yellow) are centered around the BZ center (r point) and a propeller structure consisting of the 
parabolic M -pocket (green) surrounded by four elliptical blades (magenta/brown) is centered around the BZ corner (A/-point). Top middle: 
FS topology of the 1st BZ (coloring of bands equal to left panel) together with the node of the extended s-wave gap (red). Top right: Contour 
plot of the magnitude of the extended s-wave gap given by Eq. ( |27} with Ao = 10 meV. In the middle (BZ center) the gap exhibits its 
maximum (+Ao) and in the corner it reaches its minimum (— Ao). Bottom right: Absolute value of the extended s-wave gap at the FS for 
the different bands (coloring equal to left and top middle panel) as a function of angle <f> (origin of polar coordinates always in the middle of 
the band). The angular variation of the gap on the inner V parabola (blue) and M -pocket (green) is equal to the one of the outer Y parabola 
(yellow) but not resolved on this energy scale. The elliptical blade is transformed to a parabolic band for this illustration. Color online. 




1 2 3 4 5 
Raman shift [A2] 

FIG. 4: Unscreened (black) and screened (red) Raman response in 
A\ g symmetry of the entire 1st BZ shown in Fig.|3]with two different 
constant SC gaps. The square root singularity of the smaller gap 
(present on the outer T-parabola) is removed by screening while the 
one at the bigger gap (present on all other bands) remains. Color 
online. 



C. A\ g Raman response with an extended s-wave gap 

A next step in the calculation of Raman responses for the 
chosen band structure is to drop the assumption of constant 



SC gaps and to use one general /c-dependent SC gap that de- 
termines the gap magnitude on each band dependent on its 
location in the 1st BZ. Experimental results are at this point 
still controversial even for compounds of the same stoichio- 
metric family (122 in this case) with evidence for gap nodes 
in some experiments 19 and meas urem ents of finite gaps hav- 
ing no nodes in other experiments^ 2 ^. However the candida te 
favored by the majority of publications, for example^ 2 * 21 * 22 ! is 
the extended s-wave gap which is represented by the follow- 
ing analytical expressiorPl 

A(k) = A (cos(7r/c 2: ) + cos(7rfc a )) . (27) 

Here k runs from —1 to 1. This gap exhibits a diamond shaped 
node in the 1st BZ, shown as a red line in the top middle panel 
of Fig. [3] which does not coincide with any of the bands at 
the FS. Thus all bands are fully gaped but there is a varia- 
tion of the gap magnitude on either one of them. This angular 
variation on each band at the FS (after a transformation into 
polar coordinates) is displayed in the bottom right panel of 
Fig. [5] The extended s-wave gap has its maximum Ao at the 
T-point and its minimum of the same magnitude but negative 
sign (— Ao) is lo cated at the M-point. Measured values of 
the maximum gap 6 l 23 l 24 lrange in between 6 meV and 12 meV 
supporting the choice of Ao = 10 meV for the following cal- 
culations. 

Combining Eq. ([3]l and Eq. Q at T = and noting that 
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bi, Mp, n, 



hole-hole 



hole-electron 




10 20 30 

Raman shift [meV] 



FIG. 5: Unscreened single band Raman responses in Ai g symmetry 
of all free electron bands of the band structure displayed in Fig. [3] 
(same coloring) with extended s-wave gap. All bands show a thresh- 
old at twice the minimal gap value and a log-singularity at twice the 
maximum gap value. These features are only resolved on a smaller 
energy scale for the inner T parabola (yellow) and the M -pocket 
(green) because of the smaller variation of the gap magnitude (see 
Fig. |3j. The blade exhibits an additional discontinuous falloff at 
higher energies corresponding to the top of the band in the SC state. 
Color online. 



the Raman vertex is a constant in the applied model leads to 
/ \ 1 2 4A2(k) 



in E n (k) [4£2( k )_( w + ia )2 ] 



(28) 



for the Raman susceptibility of the n th band. For the nu- 
merical calculations the sum is replaced by a 2D A; -integration 
in polar coordinates and a constant broadening parameter of 
a = 0.01 meV representing a clean system is used. The 
obtained unscreened single band Raman responses are dis- 
played in Fig. [5] in the corresponding coloring to the 1st BZ 
FS in Fig. [3] The plots show three main features that are 
seen best in the magenta response of the blade since this is 
the band exhibiting the largest gap variation at the FS. There 
is a discontinuous falloff at about 25 meV corresponding to 

2, 



A 2 

bl. center 



which denotes the top of the band in 

the SC state above which no quasi particles can be excited. 
This feature is seen at higher energies for the other bands 
(since they have larger fi's) and thus is not present in this plot. 
At twice the maximum SC gap there is a fog-singularity and 
at twice the minimum SC gap a threshold can be seen. These 
two features are also obtained for the inner T-parabola (yel- 
low) but can only be found at a much smaller energy resolu- 
tion for the other two bands since those have a much smaller 
gap variation (see lower right of Fig. [3]). The nature of the 
Zog-singularity and the threshold are found in analytical inves- 
tigations of the imaginary part of Eq. < (28] > using a simplified 
gap of the form A(</>) = A G // + Ao cos(40) with (f> running 



*£ 

u 




10 20 30 10 20 30 

4 Bl + Mp + n + 




10 20 30 

Raman shift [meV] 

FIG. 6: Screened multiband Raman responses in A\ g symmetry with 
extended s-wave gap for combinations of the free electron bands 
shown in Fig. [5] Top left: Raman responses for all combinations 
of two hole-like bands (4 Bl & To: black, 4 Bl & Ti: red, Ti & 
To: orange). Top right: Raman responses for all combinations of a 
hole-like and an electron-like band (Alp & 4 Bl: black, Mp & Ti: 
red, Mp & To: orange). Bottom: Raman response of all four bands 
of the model. The two band terms shown in the two upper panels 
renormalized by all four bands sum up to this response resulting in 
a flat continuum above twice the minimum gap of the blade and a 
sharp peak near twice the maximum gap of the M -pocket. Since 
all extremal gap values are different for all bands all singularities are 
removed for the screened two and four band responses. Color online. 



from to 27T. Noting that the integration must only run over 
the real part of the integrand leads to differences in the inte- 
gration boundaries at twice the minimum and maximum SC 
gap resulting in the two distinct features (for a detailed inves- 
tigation seeU^I). 

Knowing the single band Raman response functions is now 
used to calculate all possible combinations of the screened 
A\ g two band Raman responses and the screened A\ g Ra- 



man response of all bands. According to Eq. (19i the for- 
mer ones allow to explain the origin of the features seen in 
the latter one since it consist of a sum of the renormalized two 
band terms. In the upper left of Fig. |6]the three combinations 
of two hole-like bands are shown and in the upper right the 
corresponding plot for the combinations of the three hole-like 
bands with the electron-like M -pocket is illustrated. The total 
screened A\ g Raman response, containing all bands in the 1st 
BZ, is displayed in the bottom panel of Fig. [6] Its main fea- 
ture is a sharp peak at about 20 meV whose origin lies in the 
three electron-hole two band responses. They all show a peak 
slightly below twice the maximum SC gap of the M -pocket 
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and those add up to the observed peak in the total Raman re- 
sponse. The three hole-hole two band terms generally show 
a weaker response since the screening is more effective for 
them leading to a rather featureless broad continuum down 
to the energy of the threshold of the blade. Since all bands 
have different crossing lines with the FS the maximum SC 
gaps are different on all of them and thus all Zog-singularities 
are removed corresponding to the argumentation in chapter III 
B. The B 2g Raman response again only consists of the four 
blades which results in a Raman response entirely equal to the 
magenta single band Raman response of the blade seen in Fig. 
[5]but multiplied by 4. All other symmetry channels are again 
vanishing. 



D. Comparison with experiment 

The recent measurements of Raman responses^ on 
Ba(Fei_ x Co x ) 2 As 2 only partially resemble the obtained 
Raman responses with an extended s-wave gap. While the 
negligible contribution of the SC state to the B\ g and A 2g 
response in the experimental data is consistent with the van- 
ishing response in those channels for the used band structure 
model the A\ g and B 2g responses are quite different. The 
structure in the B 2g response of the data points towards a 
power law behavior for small energies and thus supports a SC 
gap with nodes and not the threshold observed in this work. 
Even stronger disagreement is seen in the A\ g channel where 
no sharp peak appears in the data but a rather broad dome, 
peaked at higher energies than the peak in the B 2g response, 
is found. However the Zog-singularity shape of the B 2g peak 
of the here presented Raman response is also present in the 
data. One way to tune our band structure model in order to 
get closer to the measured data is to shift the blades further 
towards the node line of the gap. To be consistent with the 
measured peak positions the effective masses of all bands, ex- 
cept the outer T-parabola, and the chemical potential of the 
A/-pocket are changed and the maximum gap value is set to 
Ao = 7 meV . All this and a larger broadening parameter 
a = 0.2meV leads to the FS topology (FS*) and the Raman 
responses illustrated in Fig. [7] The threshold in the B 2g chan- 
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FIG. 7: Left: FS* of 1st BZ of the modified free electron band struc- 
ture model together with the node of the extended s-wave gap (col- 
oring of bands as in all other figures). Right: Raman responses in 
B'2 g (green) and screened A\ g (blue) symmetry. The i?2 9 response 
is multiplied by 5. Color online. 



nel (green) is now at small energies. The total shape is close 
to the experimental one but to obtain the proper intensity the 
shown Raman response is multiplied by 5 since the spectral 
weight of the B 2g channel is quite small in the calculated re- 
sponse. This is a consequence of a lack of contribution by 
other bands than the blades to the B 2g channel and might also 
be due to some resonance effects in the experiment that are not 
taken into account in the calculations within effective mass 
approximation. The A\ g Raman response (blue) in the right 
hand side of Fig. [7] is now also quite close to the measured 
one except for the still quite obvious peak at about 13 meV. 
In general the agreement between the Raman responses of the 
data on Ba(Fei- x Co x )2As 2 and the modified free electron 
band structure is quite good which shows that the obtained 
shape of the data does not necessarily have to be explained 
with gap nodes. Furthermore it is not clear that the measured 
compound is in the clean limit still leaving the possibility of 
rather strong impurity scattering that could also explain the 
low energy tail of the data. With the introduction of possible 
scattering centers into the FeAs layer (where the important 
physics for SC is assumed to take place) by the replacement 
of Fe by Co this scenario is quite reasonable. Certainly more 
data on other compounds is needed to clarify these issues. 



V. CONCLUSION 

In this work we performed model calculations of Raman re- 
sponses for a multiband superconductor in a 2D free electron 
model. For a constant SC gap parameter analytical calcula- 
tions were used to investigate multiband effects in the screen- 
ing correction term in the fully symmetric A\ g channel. It 
was shown that for two parabolic bands with equal SC gaps 
and opposite signs of effective masses the screening term van- 
ishes which also was proven to hold true for two elliptical 
bands with equal ellipticity. Moreover the removal of the sin- 
gularities through screening for two bands with different SC 
gap sizes was presented. This feature remained in the Ra- 
man response of multiple bands for all contributing bands ex- 
hibiting a uniquely present SC gap size and the example of a 
four band response was discussed. Additionally is was shown 
through analytical calculation that multiband responses con- 
sist of a sum of two band terms normalized by all contributing 
bands. 

The general effects found for a constant SC gap parame- 
ter were then used to explain the SC Raman response of a 
model free electron band structure inspired by ARPES mea- 
surements on iron-pnictide materials. Using a SC gap of ex- 
tended s-wave symmetry the single band response features, 
a threshold at twice the minimum SC gap value and a log- 
singularity at twice the maximum SC gap value, were clari- 
fied. Furthermore we identified that the main features of the 
four band Raman response of the model band structure with 
the extended s-wave gap were originating from the two band 
terms consisting of bands with opposite signs in the effective 
masses. A comparison of this four band Raman response with 
experimental data was done and by modifying the band struc- 
ture parameters an acceptable agreement was demonstrated. 
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